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distribution is infinitely elongated along one direction, is considered here. The energy loss is 
shown to be strongly time and directionally dependent and its magnitude can much exceed the 
energy loss in equilibrium plasma. 
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1. Introduction 

When a highly energetic parton travels through a quark-gluon plasma (QGP), it losses its 
energy due to, in particular, elastic interactions with plasma constituents. This is called collisional 
energy loss and was computed for equilibrium QGP twenty years ago, see the review [1] and the 
handbook [2]. The quark-gluon plasma produced in relativistic heavy-ion collisions, however, 
reaches a state of local equilibrium only after a short but finite time interval, and during this period 
the momentum distribution of plasma partons is anisotropic. It is important to note that a plasma 
with an anisotropic momentum distribution is unstable (for a review see [4]). Collisional energy 
loss has been computed for an anisotropic QGP in Ref. [3], but the fact that unstable systems are 
explicitly time dependent as unstable modes exponentially grow in time was not taken into account. 

We have developed an approach, see [5, 6] for a preliminary account, where energy loss is 
studied as an initial value problem. The approach is applicable to plasma systems evolving quickly 
in time. We compute the energy loss by treating the parton as an energetic classical particle with 
SU(A r c ) color charge. For an equilibrium plasma the known result is recovered and for an unstable 
plasma the energy loss is shown to have contributions which exponentially grow in time. In Refs. 
[5, 6] we have calculated the energy loss in a two-stream system which is unstable due to longi- 
tudinal chromoelectric modes and found that it manifests strong time and directional dependence. 
In this paper we focus on an extremely prolate quark-gluon plasma with momentum distribution 
infinitely elongated in one (beam) direction. Such a system is unstable due to transverse chro- 
momagnetic modes and the spectrum of collective excitations can be obtained in explicit analytic 
form. The system has thus nontrivial dynamics but the computation of energy loss is relatively 
simple. After a brief presentation of our approach, we show some of our results. The energy loss 
grows exponentially and after some time its magnitude is much bigger than in equilibrium plasma. 

2. Formalism 

Our approach is classical and thus we start with the Wong equations [7] describing the motion 
of classical parton in a chromodynamic field. 

2.1 General energy-loss formula 

The Wong equations [7] read 

^(t), (2.1) 

gQ a (T)Fr(x(T))u v (r), (2.2) 

-sr%(T)<(*(T))<2 c (T), (2.3) 

where X, x^(t), u^(t) and p^{x) are, respectively, the parton's proper time, its trajectory, four- 
velocity and four-momentum; Fa and A% denote the chromodynamic field strength tensor and 
four-potential along the parton's trajectory and Q a is the classical color charge of the parton; g is 
the coupling constant and a s = g 2 /4n is assumed to be small. We also assume that the potential 
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vanishes along the parton's trajectory i.e. our gauge condition is m^(t)Aq (x(t)) = 0. Due to 
Eq. (2.3) the classical parton's charge Q c (i) is a constant of motion within the chosen gauge. 

The energy loss is given directly by Eq. (2.2) with jj. = 0. Using the time t = yr instead 
of the proper time X and replacing the strength tensor Fa by the chromoelectric E„(f,r) and 
chromomagnetic B„(?,r) fields, Eq. (2.2) gives 

dE(t) 



dt 



gQ"E a (t,r(t))-v, 



(2.4) 



where v is the parton's velocity. Since we consider a parton which is very energetic, v is assumed 
to be constant and v 2 = 1, but the parton's momentum and energy vary. 

Since we deal with an initial value problem, we apply to the field and current not the usual 
Fourier transformation but the one-sided Fourier transformation defined as 



/(<B,k) = £dt J d're^'-^fit^), 
n \ r +i ° d(0 f d3fi -i(m-k-r) w , \ 



(2.5) 
(2.6) 



where the real parameter a > is chosen is such a way that the integral over co is taken along a 
straight line in the complex ft)— plane, parallel to the real axis, above all singularities of f(co,k). 
Introducing the current generated by the parton } a {t,r) = gQ"\5^ 3 \r — \t), and using Eqs. (2.5, 
2.6), Eq. (2.4) can be rewritten: 



dE{t) 
dt 



(2.7) 



where co = k ■ v. 

The next step is to compute the chromoelectric field E a . Applying the one-sided Fourier 
transformation to the linearized Yang-Mills equations, which represent QCD in the Hard Loop 
approximation, we get the chromoelectric field given as 

£<(t»,k) = — j(E _1 )^(ca,k) [coji(co,k) + ej k, k k B l 0a (k) - ftx<(k)] , (2.8) 

where Bo and Do are the initial values of the chromomagnetic field and the chromoelectric induc- 
tion, and D' a (co,k) = e' } \(0,k)El(co,k) with £ ,y (ft),k) being chromodielectric tensor which equals 

^ 2 ' d 3 p v l df(j>) 



s' J (co,k) = d lJ + 



1 



2co J (2k) 3 ft>-k-v + /0+ dp k L\ co / co 
where /(p) is the momentum distribution of plasma constituents. The color indices a, bare dropped 
as e(ft),k) is a unit matrix in color space. The matrix £ !i (ft),k) in Eq. (2.8) is defined 

I i7 (ft),k) = -k 2 8 ij + k i k j + co 2 e ij {co,k). (2.9) 

Substituting the expression (2.8) into Eq. (2.7), we get the formula 



dE(t) 
dt 



8Q V j 
icogQ"v ] 



dco 

oo-\-i(j 2,7li 



+ £ J kl k% a (k)-coDi a (k) 



(2.10) 



CO- CO 

As seen, the integral over co is controlled by the poles of the matrix (ft),k) which represent the 
collective modes of the system. 
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Figure 1: The parton energy loss per unit time in equilibrium plasma as a function of 

2.2 Equilibrium plasma 

When the plasma is stable, all modes are damped and the poles of £~ 1 (ft), k) are located in the 
lower half-plane of complex CO. Consequently, the contributions to the energy loss corresponding 
to the poles of (ft),k) exponentially decay in time. The only stationary contribution is given by 
the pole co = CO = k • v. Therefore, the terms in Eq. (2. 10), which depend on the initial values of 
the fields, are neglected and Eq. (2.10) provides 



dE 
~dt 



-ig 2 CR 



d 3 k CO 



(2tt) 3 k 2 



1 



2, 7 2 



+ 



k 2 v 



co z 



£ L (co,k) co 2 £ T (co,k) -k 2 



(2.11) 



N 1 

where the color factor Cr, which equals for a quark and N c for a gluon, results from the 
averaging over colors of the test parton. The formula (2.11) agrees with the standard energy loss 
due to soft collisions in equilibrium QGP [2]. 

To compare the energy loss in an unstable plasma to that in an equilibrium one, we have 
computed the integral in Eq. (2.11) numerically using cylindrical coordinates, which will also be 
used for the prolate system. Since the integral is known to be logarithmically divergent, it has been 
taken over a finite domain such that — k max < ki < fc max and < kj < k max . The energy loss in an 
equilibrium plasma of massless constituents can be expressed through the Debye mass, which we 
write as 

d 3 p /(p) 



M 2 =, 



(2«y- 



(2.12) 



We define dimensionless variables by scaling dimensionful quantities by the Debye mass. In Fig. 1 
we show the energy loss in equilibrium QGP divided by g 2 ^ 2 as a function of computed for 
g = landC R = N c = 3. 



2.3 Unstable plasma 

When the plasma is unstable, the matrix (co,k) has poles in the upper half-plane of com- 
plex co, and the contributions to the energy loss from these poles grow exponentially in time. The 
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terms in Eq. (2. 10) which depend on the initial values of the fields D and B cannot be neglected, as 
they exponentially grow in time. Using the linearized Yang-Mills equations, the initial values Bo 
and Do are expressed through the current and we obtain 



dE(t) 
dt 



g 2 C R Vv l 



-oo+jct 271 J (2k) 3 



(2.13) 



CO 



8» 



(k j k k -k 2 8 jk )(^- l ) k '(m,k) + come J \a),k)(Z- l ) k \a),k) 



(0 — (O 

which gives the energy loss of a parton in an unstable quark-gluon plasma. 
2.4 Inversion of £ 

When the anisotropy of the momentum distribution of plasma constituents is controlled by a 
single (unit) vector n, it is not difficult to invert the matrix L. Following [9] , we introduce the vector 
nr defined as 

k'kK 



ri T = (S ij 

and we use the basis of four symmetric tensors 



(2.14) 



A' J ( k ) = §' J _ _ , b-j ( k ) = _ ? & (k, n) = 



D lJ (k,n) = k'n J T +k J n l T . (2.15) 



Since the matrix £ is symmetric, it can be decomposed as £ = a A + bB + cC + dD, where the 
coefficients a, b, c and d are found from the equations 

k l l} j k j = k 2 b, rf T l} j n{ = n 2 T (a + c), n^L ij k j = n|k 2 J, TrI = 2a + b + c. (2.16) 

The inverse matrix is found to be 



. _ 1 -a(a + c)B + (-d 2 k 2 nl+bc)C + adD 
2-i — A ~\~ 



a(d 2 k 2 n 2 -b{a + c)) 

and consequently, the poles of the matrix are given by the following dispersion equations 



(2.17) 



a = 0, d 2 k 2 n 2 T -b(a + c) = 0. 



(2.18) 



Writing down the energy-loss formula (2.13) in terms of the projectors A,B,C,D we obtain 
the form 



dE(t) , 



°+ ia dco f d 3 k 



o+ic 2%i J (2n) 3 



,-i(a)-w)t 



1 



-A + 



-a(a + c)B + (-d 2 k 2 nf + bc)C + adD\ 



a(d 2 k 2 n 2 -b(a + c)) 



CO 



CO (0 + (0 ^-( l A (-^ 2n T+ b ^) C + ddAD \ 

co-co co a> \a d(d 2 k 2 nj. -b(a + c)) ) 



(2.19) 



where the coefficients a,b,c,d are computed at (ft),k) and the coefficients a,b,c,d at (co,k). The 
formula (2.19) will be used in the subsequent section to compute the energy loss in an extremely 
prolate QGP. 
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3. Extremely prolate plasma 

Anisotropy is a generic feature of the parton momentum distribution in heavy-ion collisions. 
At the early stage, when partons emerge from the incoming nucleons, the momentum distribution is 
strongly elongated along the beam - it has a prolate shape with the average transverse momentum 
being much smaller than the average longitudinal one. Due to free streaming, the distribution 
evolves in the local rest frame to a form which is squeezed along the beam - it has oblate shape 
with the average transverse momentum being much larger than the average longitudinal one. We 
consider here the extremely prolate momentum distribution which is either infinitely elongated in 
the longitudinal direction (defined be the unit vector n) or infinitely squeezed in the transverse 
directions. In this case the spectrum of collective excitations can be found analytically. 



3.1 Momentum distribution 

The extremely prolate momentum distribution can be written as 

/(p) = S(p 2 -(p-n) 2 )%-n), (3.1) 

where h(x) is any positive even function such that / d 3 f(p) is finite. The integral in Eq. (2.12) 
can be used to define a mass parameter for either isotropic or anisotropic momentum distributions, 
although in the later case pt~ { cannot be interpreted as the screening length. Applying the prolate 
distribution a different mass parameter denoted as m naturally appears. It is related to pt as m = 

Since the velocity v of a massless parton as given by the distribution (3. 1) is v = n for p • n > 
and v = — n for p • n < 0, the matrix £ defined by Eq. (2.9) is found to be 



Z' 7 (co, k) = (co 2 - m 2 - k 2 )8 ij + k'k j _ W ~ ' N7 (JtV + nV) (3.2) 

'"-1 ^ — it ,nU 



m 2 k • n 
co 2 - (k-n) : 

m 2 (co 2 + (k-n) 2 )(k 2 -co 2 ) 
(co 2 - (k-n) 2 ) 2 



and the coefficients a,b,c,d are 



a(co,k) = co — m — k , (3.3) 

M co,k) = co 2 -„ 2 - 2 ft U) l - m2K + (k - n)2)(k2 2 "" 2) ^, (3-4) 

co 2 -(k-n) 2 (co 2 -(k-n) 2 ) 2 k 2 

m 2 (co 2 + (k-n) 2 )(£ 2 -co 2 ) /(k-n) 2 \ 

c(W ' k) = (co 2 -(k-n) 2 ) 2 V ' (3 " 5) 

- ~ ft^ -- 2 ^ + ( k -) 2 )( k2 -^(^, (3.6) 
fi) 2 -(k-n) 2 (co 2 - (k-n) 2 ) 2 k 2 
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1 dE 




Figure 2: The parton energy loss per unit time as a function of time for three angles © between the parton's 
velocity v and the axis z. The red points correspond to © = 0, the blue ones to = n/12 and the purple 
points to © — k/6. The solid lines represent the exponential fits to the computed points. 



3.2 Collective excitations 

The dispersion relations are obtained from Eqs. (2. 1 8) with the coefficients given by Eqs. (3.3- 
3.6). The first of these equations provides ft) 2 (k) = m 2 + k 2 . Although the second equation looks 
rather complicated, it has three relatively simple solutions 



ft) 2 (k) = m 2 + (k-n) 2 , 



(3.7) 



o 2 (k) = ^k 2 + (k-n) 2 ±^/k 4 + (k-n) 4 + 4m 2 k 2 -4m 2 (k-n) 2 -2k 2 (k-n) 2 ) , (3. 



8) 



which hold under the condition ft) 2 / (k • n) 2 . The modes (0\, 00i and ft) + , are always stable. The 
solution CO 2 is negative when m 2 k 2 > m 2 (k • n) 2 + k 2 (k • n) 2 . Writing ft) 2 = — y 2 , < y G R, 
the solutions are =h'y. The first is the Weibel unstable mode and the second is an overdamped 
mode. Denoting the angle between k and n as 6, the condition for the existence of an instability 
is cos 2 6 < m ™ +k 2 ■ One can show that for fixed k 2 the instability growth y is maximal when kin. 
Collective excitations in the extremely prolate QGP were earlier studied in [10] using a method 
different than ours. 

Using Eqs. (3.7) and (3.8) one can write 



a{d 2 \i 2 n 2 T -b(a + c)) 



ft) 2 (ft) 2 - co 2 (k)) (co 2 - of (k)) (CO 2 - (O 2 (k)) (CO 2 - CO 2 (k)) 



co- 



(k-n) 



(3.9) 



and from Eq. (2.19) one therefore sees that the energy loss in the extremely prolate system is 
controlled by the double pole at ft) = and 8 single poles: ft) = ±C0i, ft) = ±0>2, CO = ±ft) + , 
ft) = ±ft)_. Since the collective modes are known analytically, the integral over ft) in Eq. (2.19) 
can be computed analytically as well. The remaining integrals are performed numerically. 
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3.3 Numerical results 

To compute the integral over k in Eq. (2.19), we use cylindrical coordinates with the z axis 
along the vector n. Since the integral is divergent (as is the case in equilibrium (2.1 1)), we choose 
a finite domain such that —k max <k^< k max and < kj < k max with k max = 5^1. The values of 
remaining parameters are: g = 1, Cr = N c = 3. In Fig. 2 we show the parton's energy loss per unit 
time as a function of time for three different orientations of the parton's velocity v with respect to 
the z axis. The energy loss manifests a strong directional dependence and it exponentially grows in 
time, which indicates the effect of the unstable modes. After a sufficiently long time, the magnitude 
of energy loss much exceeds that in equilibrium plasma which, as shown in Fig. 1, equals 0.18 for 

^•max — 5/X. 

4. Conclusions 

We have developed a formalism where the energy loss of a fast parton in a plasma medium 
is found as the solution of an initial value problem. The formalism, which allows one to obtain 
the energy loss in an unstable plasma, is applied to an extremely prolate quark-gluon plasma with 
momentum distribution infinitely elongated in the z direction. This system is unstable due to chro- 
momagnetic transverse modes. The energy loss per unit length of a highly energetic parton is not 
a constant, as in the equilibrium plasma, but it exponentially grows in time and exhibits a strong 
directional dependence. The magnitude of the energy loss can much exceed the equilibrium value. 
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